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We report on the progress in the computation of the /3-functions of c/> 4 theory and 
QCD in the large N expansion. For the former we give an analytic formula for 
the critical exponent which encodes higher order coefficients in the series beyond 
those currently known in the MS scheme at 0(1/N 2 ) in a large N expansion. This 
allows us to deduce new coefficients in the field anomalous dimension which relate 
to Kreimer's knot theory of the divergences of a quantum field theory. A similar 
exponent is computed for the /3-function of QCD at 0(1/ Nj). 



1 Introduction 

We report on two areas where progress has been made into revealing the higher 
order structure of /3-functions using the large TV expansion. First, we recall 
that a highlight of the Pisa workshop was the breakthrough by Broadhurst in 
computing those diagrams of the six and seven loop /3-function in </> theory 
which are free of subdivergencescl This was in response to Kreimer's association 
of the divergences of multiloop Feynman diagrams involving transcendental 
numbers with knot theory and in particular positive braidsB For example, the 
appearance of the Ricmann zeta series, £(2n — 3) n > 3, is well known in higi 
order MS calculations and these were associated with (2n — 3, 2) torus knots 
The six and seven loop terms of the /3-function contained new transcendentals 
These were represented by double or triple sums which could not be reduced 
to products of C( n )- They were associated with the torus knot IO124 and the 
hyperbolic knots IO139, IO152 and H353. As these interesting knot numbers 
appear at very high orders it is important to have an efficient way of accessing 
them. The large N approach is a technique which facilitates thisla For instance 
in </> 4 theory with an O(N) symmetry the field anomalous dimension, 7(5), has 
been computed analytically in d-dimensions at 0(1/ N 3 ) at a fixed point, g c , 
using a conformal bootstrap methods The result involves a particular two loop 
Feynman integral whose e-expansion, where d = 4 — 2e, contains irreducible 
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double sums. These therefore can be associated with the more complicated 
positive braids. To find at what order in perturbation theory they arise and to 
determine the explicit MS coefficients in 7(5) associated with them one needs 
to solve the critical renormalization group relation r\ — 2 r y{g c ). Necessary for 
this is the location of the critical coupling g c = g c (e, N) which is dimensionless 
in 4-dimensions. In this article we report the result of the 0(1/N 2 ) calculation 
of f3'(g c ) and therefore g c which allows us to proceed. The remainder of the 
paper is devoted to a similar calculation of the /3-function of QCD0 though at 
0(1/Nf), which will provide new MS coefficients and therefore an independent 
check for future 4 and higher loop calculations, albeit at the easier end of the 
range of perturbative calculability. 

2 tp 4 /3-function 

The large N method to compute all orders coefficients of the renormaliza- 
tion group-iunctions at a fixed point of the /3-function has been discussed 
previouslyBQ For </> 4 theory with an O(N) symmetry we use the lagrangian in 
the form 

L = IW 1 ) 2 + - a 2 /2g (1) 

where a is an auxiliary field. This form is used so that uniqueness can be 
applied and-also to emphasis the equivalence with the O(N) a model in (2 — 2e) 
dimensionsu This <i-dimensional equivalence means that the critical exponents 
computed in either model are the same. Although various a model exponents 
aside from r\ are known at O(1/7V 2 )0 the exponent to which we define as ui = — 
(3'{g c )/2 has not been determined. It plays a role in understanding corrections 
to scaling and as it is not related through scaling laws it must be computed 
independently. 

In the large N formalisrj] this is achieved by considering corrections to 
asymptotic scaling and then solving the tf> and a Schwinger Dyson equations 
self consistently. At leading order due to the structure of the A^-dependence 
of the 2x2 matrix yielding u) there is one 2 loop and one 3 loop Feynman 
diagram to compute^ This situation is parallel to the same computation in the 
Gross Neveu modelP For the 0(1/N 2 ) calculation, in addition to computing 
the 0(1/N) corrections to the previous graphs, there are two 2-loop, eight 
3-loop, eleven 4-loup, ten 5-loop and two 6-loop integrals which need to be 
evaluated. We findpl 

wi = (2 M -1)V (2) 
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where r?i = - 4r(2/i-2)/[r(2-/r(/i-l)r(/i-2)r(/i+l)],_<J<(/ = ip(2fj,-3) 
+ ^(3 - / - - 1) - V(l), e(/ = /(/i - 1) - /(l), $(/ = #(2/* - 3) 
- V'(3 - /i) - ^'(/i - 1) + /(I) and uj = YsZo^i/N 1 with ^0 = ^-2. The 
e-expansion of Eq. (3) near d = 4 — 2e dimensions gives total agreement with 
the coefficients of the explicit 5-loop MS resultsli3 So writing the /3-function in 
the form where only the 0(1/N) and 0(l/N 2 ) coefficients are included, 



/3(g) = \{d-A)g + (a 1 N + b 1 )g 2 + Y(a r N + b r )N r - 2 g r+1 



(4) 



we find the new MS coefficients 

a 6 = [29 + 528C(3) -432C(4)]/1866240 

a 7 = [61 + 80C(3) + 1584((4) - 1728C(5)]/26873856 

b 6 = [ - 28160C(7) + 95200C(6) - 150336C(5) - 6912C(4)C(3) 
+ 14112C(4) + 24064/(3) + 11880C(3) + 5661]/466560 

b 7 = [ - 7795200C(8) + 32724480C(7) - 43286400C(6) - 3993600C(5)C(3) 
+ 31998720C(5) - 622080C 2 (4) + 8663040C(4)C(3) + 2538432C(4) 
- 11381760C 2 (3) - 7461168C(3) - 1125439]/403107840 (5) 

Knowledge of Eq. (3) means that g c is available at 0(1/ N 3 ) to all orders in 
e and therefore we can proceed in deducing coefficients in the series of 7(5), in 
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a form similar to Eq. (4). The expression for r\ — 2j(g c ) at 0(l/iV 3 ), however, 
involves the e-expansion of a 2-loop self energy integral whose evaluation in 
terms of standard integration techniques for massless Feynman diagrams such 
as uniqueness is not possible for arbitrary da Using group theory, for exarrL: 
pie, it has recently been written in terms of an 3-F2 hypergeometric functional 
Therefore, writing the perturbative form of 7(g) as 

oc 

7(5) = Y.^ N2+d <- N + e r) N ^ 2 9 r+1 (6) 
r=l 

with ei = where we have only indicated the leading three orders in the 1/N 
expansion, we deduce that the first appearance of a non-zeta number at this 
order in 1/N is in eg, or ninth loop. Explicitly 

e 8 = [1044119552C(9) - 2854457088C(8) + 2246396928C(7) 

+ 264929280C(6)C(3) + 572052480C(6) + 77856768C(5)C(4) 
- 663797760C(5)C(3) - 887356416C(5) - 127844352C 2 (4) 
+ 309288960C(4)C(3) - 476930304C(4) - 8978432<C 3 (3) 
+ 286906368C 2 (3) + 677644800C(3) 

+ 12386304£/ 62 - 96663009] /216710774784 (7) 

In Eq. (7) Uq2 is the level 8 non-zeta number first discovered in Ref. [11]. In 
the /3-function calculation it occurs at sixth loop and is associated with the 
819 = (4, 3) torus knotS 

3 QCD /3-function 

We now report on a similar calculatioii for the /3-function of QCD but at 
0(1/Nf) where Nf is the number of quark flavours. An essential feature of the 
previous section was the equivalence of 4> theory at the fixed point to another 
model. The property of universality meant that exponents calculated in cither 
model gave the same results. This principle can be exploited in the case of 
QCD. In large Nf it has been demonstrated, at_least in leading order, that 
QCD is equivalent to a lower dimensional modelEj Therefore it is possible to 
compute d-dimensional critical exponents in either model but the coefficients 
of the perturbative renormalization group functions of QCD are deduced from 
knowledge of the location of the QCD fixed point. The main advantage of 
using this equivalence in this instance is that the lower dimensional model 
has a simpler form. In particular it has fewer interactions which immediately 
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reduces the number of Fcynman diagrams which need to be computed. More 
concretely the lagrangian for QCD, in standard notation, is 



L = i^W 1 ' %^ (8) 

The model to which it is equivalent as Nf — > oo is the non-abelian Thirring 
model, (NATM) , which is perturbatively renormalizable in strictly two dimen- 
sions and has the lagrangian, 

(A a ) 2 

L = ifi 1 ^ 1 - K -*L- (9) 

where the coupling constant A is dimensionless in 2-dimensions. Clearly the 
NATM contains only a single quark gluon interaction. In using Eq. (9) to 
compute exponents there would appear to be no contributions to quantities in 
QCD which are built with the triple and quartic gluon self-interactions. It was 
demonstrated, however^ that these pieces are contained in subdiagrams which 
involve a quark loop which connects three and four gluon lines respectively. 
In other words using the NATM with the critical point massless propagators 
diagrams with such loops reproduce the contribution to the perturbation the- 
ory in QCD of both the loop itself and also the gluon self-interactions. Indeed 
this is apparent in the results we have previously computed using Eq. (9) in 
relation to the known explicit 3-loop MS resultsQ Further in dealing with non- 
abelian gauge theories we use a covariant gauge fixing. Therefore ghost fields 
must also be included in both lagrangians. It turns out that since a loop of 
ghost fields gives one power less of Nf than at leading order in 1/A^ they do 
not contribute to the exponents. 

To find to one computes the dimension of the operator associated with that 
coupling^ From Eq. (8) this is G = (G°„) 2 and from the last term of Eq. (8) it 
is easy to deduce the scaling law which relates u> to the anomalous dimension of 
G and the gluon anomalous dimension . As the 0{\/Nf) result for the abelian 
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sector has been reported previously 3 ! 1 all that needs to be computed are the 



diagrams in the gluon 2-point function with G inserted whose colour factor is 
proportional to C2(G). There are two such graphs; one two loop and one three 
loop. Calculating these with the massless critical point propagators we obtain 
the result 



(2/i-3)(/x-3)C 2 ( J R) 

(4ju 4 - 18/i 3 + 44/x 2 - 45^ + 14)C 2 (G) 
4(2/x-l)( M -l) 



o 



m 



T(R) 



(10) 



5 



This has been computed in an arbitrary covariant gauge and we have checked 
that the gauge parameter vanishes as it ought. Except in the Landau gauge, 
however, there is potential for the operator G to mix with (d^A") 2 . It turns 
out that the mixing matrix at criticality is in fact diagonal. The final check we 
have is in relation to the e-expansion..aLEa. (10) and we record that it agrees 
with the 3-loop MS QCD /3-functiont3'EjO This agreement further justifies our 
previous remarks concerning the use of the NATM in computing exponents in 
QCD. In similar notation to Eq. (4) we deduce the new higher order coefficients 

a 4 = - [154G 2 (i?) + 53G 2 (G)]/3888 

a 5 = [(288C(3) + 214)C 2 (i?) + (480C(3)- 229)C 2 (G)]/31104 
a 6 = [(864C(4) - 1056C(3) + 502)G 2 (i?) 

+ (1440C(4) - 1264C(3) - 453)G 2 (G)]/233280 
a 7 = [(3456C(5) - 3168C(4) - 2464C(3) + 1206)C 2 (i?) 

+ (5760C(5) - 3792C(4) - 848C(3) - 885)C 2 (G)]/1679616 (11) 

4 Conclusion 

We have provided the two leading MS numbers in the quintic in N which 
appears as the coefficient in the /3-function of 4 theory at six loops. Although 
information on the hardest diagrams has been determined in the N — 1 caseJIrB 
computation of the remaining numbers to complete the six loop term is still 
far off. For this to be achieved a very high degree of ingenuity will be required 
to produce an algorithm which can then be implemented in some symbolic 
manipulation programme. 
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